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The fundamental theorem of Riemannian geometry is inverted for analytic Christoffel symbols. 
The inversion formula, henceforth dubbed Ricardo's formula, is obtained without ancillary assump- 
tions. Even though Ricardo's formula can mathematically give the full answer, it is argued that the 
solution should be taken only up to a constant conformal factor. A procedure to obtain the Christof- 
fel symbols out of unparameterized geodesies is sketched. Thus, a complete framework to obtain 
the metric out of measurements is presented. The framework is suitable for analysis of experiments 
testing the geometrical nature of gravity. 
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I. INTRODUCTION 

The archetypical experiment in General Relativity con- 
sists in following the trajectory of tests particles and de- 
termining the presence of the gravitational field from the 
mutual deviation of these world lines. Postulating that 
test-particle trajectories are geodesies is not a trivial step. 
After all, a world line is deemed to be a geodesic only 
if its length is a extremum. The notion of length re- 
quires a metric. Therefore, the analysis of any exper- 
iment designed to test the strong equivalence principle 
should avoid making assumptions that are tantamount 
to having a metric available, even an unknown one. On 
the other hand, if the very existence of the metric is not in 
question, then there are methods to figure either the con- 
nection or the curvature out of the geodesies [l|, [3, H, |3] ■ 
From there, one can apply any of the formal methods 
available in the literature to find the metric [1, 0] . 

Physicists favor computing the curvature because the 
connection, being a collection of tensors but not a tensor 
itself, is usually not deemed to encode a physical object. 
Besides, if one had the connection, the curvature is just a 
few calculations away. Thus, finding the metric out of the 
connection or out of the curvature are equivalent steps in 
the search of the metric out of observations. Mathemat- 
ically speaking however, the connection is at the center 
of the Fundamental Theorem of Riemannian Geometry: 
For a given metric, there is one and only one torsion-free 
connection that parallel transports the metric. In a co- 
ordinate system, the components of such connection is 
given by Christoffel's formula: 

r"/37 = ^5"^ (5m/3,7 + 5m7:/3 - 5/37,m) • (1) 

The Inverse of the Fundamental Theorem of Rieman- 
nian Geometry is the solution of the above partial dif- 
ferential equation in the metric when the Christoffel 
symbols and appropriate boundary conditions are given. 
Hereafter, we will shorten this as "finding the metric" . 
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This may include the calculation of the curvature as an 
intermediate step. With "finding the Christoffel sym- 
bols" , we will refer to the step of computing the compo- 
nents of the connection out of trajectories in the coordi- 
nate system in which the trajectories are specified. 

Breaking the calculations in these two steps is implicit 
throughout the literature. Perhaps the only counterex- 
ample is SingatuUin's calculations of the metric out of 
properties of geodesies valid for static metrics [7] . 

Finding the metric is tackled informally by postulating 
the most general metric compatible with some simplify- 
ing symmetries {e.g. spherical symmetry), computing the 
connection from such metric, and best-fitting it against 
the measured values. The main drawback of this pro- 
cedure is that the existence of the metric is assumed a 
priori. These analyses should not be used to test the 
strong equivalence principle for example. 

Known to the author, there are two classes of pro- 
cedures to find the metric. The algorithm devised by 
Schmidt Q requires computing the holonomy group (32| 
and all of its quadratic invariant forms. However as 
pointed out by Ihrig, this is rather only a formal solution 
on the account that it still involves solving the parallel- 
transport equations [8]. Moreover, the emphasis on the 
holonomy group makes the problem unnecessarily diffi- 
cult because the solutions are sought independently for 
each holonomy class {e.g., as done by Hall [9.] or Ghanam 
[l0|). Also, classifying the connections in clear-cut classes 
is incompatible with the fuzzy values associated with un- 
certainties. 

The second algorithm to find the metric was originally 
presented by Ihrig ^Qj . It determines the metric using the 
curvature instead of the connection [ll|, [l^, [H, [i3l and 
the solution is found with some restrictions (for example, 
when certain matrix is invertible, c./. Eq. (14) in [13j). 
It is an attractive method in that there are no differential 
equations to solve. Still, this approach is not desirable to 
compute the metric uncertainty, because the restrictions 
are not guaranteed to be satisfied for all values inside 
of the interval of confidence (of the curvature) . None of 
these two programs deal satisfactorily with the issue of 
boundary conditions. 

This paper presents a direct method to find the met- 
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ric without assuming any simplifying symmetries. It is 
only required that the Christoffel symbols are analytic in 
a geodesically connected patch [33] and that the set of 
isometries are known at a point of such patch. The an- 
alyticity requirement is stronger than the usual smooth- 
ness condition but it should not pose practical re- 
strictions for "real-life" experiments. The solution works 
for spaces of any dimension and signature. 

For the sake of completeness, we will review the prob- 
lem of finding the Christoffel symbols in the next sec- 
tion. This will allow us to cover the whole procedure 
from measurements to metric with the same philosophy 
of avoiding extra assumptions. The main result is de- 
rived in the third section and analyzed in the following 
one. An example is provided in the appendix. 



This paper is concerned with the general experimental 
situation where each trajectory is characterized only by 
the coordinates of its events. That is, the experimentalist 
provides us with the unparameterized trajectories of test 
particles. Since the afhne parameter is unknown, the 
Christoffel symbols cannot be read off directly from the 
geodesic equation 



+ r"^^ i^i" = 0. 



(2) 



Instead, one can parameterize the trajectories in terms 
of a coordinate, say x^. Eq. ^ is then transformed using 



x"" = x^fti^, and 

■•a a • b • b , a --b 



(3) 
(4) 



II. THE CHRISTOFFEL SYMBOLS 

The problem of finding the Christoffel symbols (or the 
curvature) out of trajectories of test objects has been 
studied before. Synge developed the five-point method 
[Hi 13 1 Pirani followed the trajectory of a cloud 0, 
Szekeres proposed the gravitational compass [3] , Ciufolini 
used Fermi- Walker transport In all of these works, 
the trajectory is known as a function of the affine pa- 
rameter. Otherwise in those derivations, it would have 
been impossible to assert that some 4- vector is the (nor- 
malized) 4-velocity of the test particle. Moreover, as 
Quevedo pointed out in relation to the above cited meth- 
ods [3l , "the measurement of the curvature is associated 
with an orthonormal tetrad comoving with an observer" . 
However, this implies an a priori notion of orthonoma- 
lity. In other words, one cannot assert that a tetrad is 
orthonormal until the metric has already been found. A 
similar case is made against using Fermi- Walker trans- 
port whose defining equation implicitly involves the met- 
ric. 

The affine parameter might be available under some 
special cases and further information can be extracted 
However, it is important to obtain the connection 
without assuming that the connection is a metric connec- 
tion because alternative theories of gravity might differ- 
entiate [13 between affine geodesies and metric geodesies 
[sl ]. In fact, deahng with geodesies to find the correct 
metric might not be trivial and errors might go unnoticed 
for generations p^ . 

If the "metricity" of the connection is not chal- 
lenged, then the uniqueness of the Christoffel symbols has 
been proven for vacuum [16], for Friedmann-Lemaitre- 
Robertson- Walker 19], and for certain Einstein 
spacetimes. As shown by Marmo et al, the requirement 
of being a metric connection is a very strong one: There 
are connections that can be derived from a Lagrangian 
but fail to be metric connections 21| . An important case 
of nonmetric connection is Newtonian mechanics pT| . 
Thus, determining the "metricity" of a connection be- 
comes a test of General Relativity. 



where x"' = x'^{x'') is the a-th coordinate considered as a 
function of the 6-th coordinate, indicates the deriva- 
tive of respect to x^, and ^j, stands for the second 
derivative. There is no summation on repeated Latin 
indices. We have then 



x\ = 1 and x\ = 



(5) 



and the geodesic equations for a and b become 



i' + rVA<.i'i' = o, (6) 

a • b • b , a --b , -pa u • b • b r, /y\ 

X j^fj XX -r X j^X -\- i ^lyX ^ X X X — U. ylj 

Assuming that the geodesic does not correspond to a par- 
ticle standing still and that x'' is not a constant of the 
motion, we can remove x^ from the equations above and 
arrive to 



,bb ■ 



n I 



b 



0, 



(8) 



which is a differential equation with x'^ as the indepen- 
dent variable and no reference to the affine parameter 
whatsoever. 

Because of the parameterization x"' = a;°(a;'') on the 
trajectories, Eq. ([5]) can be interpreted as a linear equa- 
tion on the unknown Christoffel symbols. In four di- 
mensions, each trajectory can potentially give 12 scalar 
equations by taking different values of a and b in Eq. ([5]). 
Thus, one needs to consider at least four trajectories per 
point in order to find the 40 values of F^^^^. It can 
be proven however that Eq. ([5]) does not determine the 
Christoffel symbols uniquely even if we used more than 
four trajectories trough the same point. 

The non-uniqueness of the mapping from trajectories 
to connections does not come out as a surprise. After all, 
we have Levi-Civita's theorem 123 : two connections F" 



and F ^ have the same unparameterized geodesies if and 



only if their difference can be written as 6"^ 
where is an arbitrary covariant vector. 
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III. THE METRIC 



From Eq. ([T]), it is trivial to obtain a linear first-order 
partial differential equation on the metric components: 



where 



Q/37 •— Q7" 13 +" a'^ I3j- 



(9) 



This equation can be easily solved by converting it into an 
ordinary differential equation. We can accomplish this by 
inverting the reasoning that defines a partial derivative 
as a special case of the directional derivative. We will 
solve Eq. ^ along a path 



(10) 



from the point xq (where the isometrics are known) to the 
point a; / (where the metric is to be evaluated) . Since the 
ChristofFel symbols can now be interpreted as functions 
of the parameter A, then Eq. after multiplying both 
sides by -gj, becomes an ordinary differential equation: 



dA 



(11) 



It can be solved, say using Frobenius method, and eval- 
uated at A = 1. 

The expansions of 5 := gap and T := ^^'^ ap^i^f^ ~ 
xq^) around a;o can be written as 

oc 

aiX' , and 



i=0 



i=0 



(12) 
(13) 



The tensor indices have been removed to avoid clutter in 
the expressions. After replacing Eqs. ([T^ and into 
Eq. (fTTj) and gathering the powers of A, we find 



■E 

i=0 



5 = ^Oi, 



(14a) 
(14b) 

(14c) 



i=0 



where go is the metric at a;o s-nd the operation is defined 
by p^g : 
formula. 



by p*q := Pp.u'l^'^ ap- ^^^^ dl) shall be named Ricardo's 



IV. ANALYSIS 



Finding the metric has been studied extensive!' 

, S mm d, [H, Q M, M, M, El, m, [13, 1251 



1271 . |28|. These works provide a mixture of somewhat 
general algorithms and studies of existence and unique- 
ness. Ricardo's formula has the advantage of being quite 
universal. The only requirements are analytic Christof- 
fel symbols in a geodesically connected patch. The ana- 
lyticity prerequisite excludes pathological examples {e.g. 
Eqs.(6)-(8) of [|). The demand of geodesical connected- 
ness will take care of ambiguities arising from singulari- 
ties. 

While the convergence of Eqs. (HH) can be proven for 
any set of in?(n + l)/2 analytic functions, where n is 
the dimensionality of the manifold, it has already been 
proven that not any such set of functions corresponds to 
the Christoffel symbols from a metric (an example can be 
found in Eq. (1) of [25]). Edgar [23] gave the necessary 
and sufhcient conditions for a connection to be locally 
metric [391: 



"-p(Q^ P)^S 



0, 



(15) 



where h is the so-called metric candidate. In our case, 
h is just g of Eq. (|14cp . Atkins extended Edgar's results 
when he proved that an analytic locally metric connec- 
tion on an analytic simply-connected manifold is met- 
ric ^2^. Thus, the integrability condition of Eq. ^ is 
Eq. (|15p . In practice, it is easier to just check that the 
connection is recovered using Eq. p^. 

Now, let us focus on the uniqueness of the solution. Be- 
ing a linear first-order partial differential equation with 
analytic coefhcients, the Cauchy-Kowalevski theorem im- 
plies that Eq. ^ should have a unique solution deter- 
mined by the appropriate boundary conditions (see page 
348 in [301 ). However, the solution of Eq. pTjl might be 
path dependent. We can prove that this is not the case 
by integrating Eq. pT|) along a closed loop. The result- 
ing right-hand side can be cast as the surface integral of 
f/i(a'^''/3)7(5' which is known to vanish (the connectedness 
requirement avoids enclosing any singularity within the 
integration path). Therefore, J)' dA = and the metric 
defined by Eqs. (|14p is independent of the path chosen in 

Eq. uni). 

For physicists, it is reassuring to know that the signa- 
ture of the metric does not change from event to event in 
spacetime. In order to prove this, it suffices to compare 
the signature of the metric at Xf and the signature of go. 
Since the metric given by Eq. (|12p is continuous, the sig- 
nature at Xq can be different from the signature at only 
if there exists at least one point in a path joining xq and 
Xf where the determinant of the metric vanishes. But at 
such point, because of Eq. ([1]), the Christoffel symbols 
would cease to be analytic. Note that the constancy of 
the signature is a consequence of the analyticity require- 
ment and it is valid for geodesically connected manifolds. 
In particular, it has nothing to do with the equations of 
motion that yielded the geodesies. 

The procedure in the previous section has broader ap- 
plication to solve a larger class of partial differential equa- 
tions. Also, there are implications in the handling of the 
boundary conditions. Both of these issues are out of the 
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scope of this paper. For the purposes of this article, we 
note that the boundary conditions have been replaced by 
the knowledge of the set of isometrics at xq . This means 
that we know (^g, the actual value of the metric at xq, 
used in Eq. (|14a[) up to a positive factor. This factor 
propagates to all terms of the series of Eq. (I14cp because 
of Eq. (|14b[) . Therefore, our solution is given up to a 
positive constant conformal factor. If, instead of using 
the set of isometrics as boundary conditions, we used 50 
then the conformal factor would have been fixed. 

In the literature [l6l.[26j. the indeterminacy of the over- 
all factor in the metric is usually taken to correspond 
to the freedom of choosing the unit of length. This is 
only partially correct. The coordinates of the events in 
the geodesies were determined with some measurement 
apparatus. The calibration of the apparatus settles the 
units. 

What is missing is a way to relate the units for the 
various axis. For example, why do we stipulate that a 
light year along the x direction is the same as a the light 
year along the y direction or the same as a year along 
the t direction? Such relationships need to be specified 
only at one point. The connection then translates the 
relationships to the rest of the patch. The relationships 
are statements about the measurement apparatus and 
they are independent of the Physics that determined the 
geodesies. For instance. Local Lorentz Invariance cannot 
be used to lift the indeterminacy. 

The relationships gather in the set of isometrics but 
not on the abstract group of isometrics. In order to un- 
derstand why we focus on the set of isometrics, consider 
the role of c on 

ds^ = -c^di^ + dx^. (16) 

Different values of c correspond to the speed of light be- 
ing measured using different units. The set of isometrics 
is not the same for each value of c (in the same sense that 
the interval [0,1] is different from the interval [0,2]), al- 
though the abstract group is the same. Pointwise, the 
value of the metric determines the set of isometrics but 
the set of isometrics decides the metric only up to a fac- 
tor. In short, the measurement apparatus fixes the units 
and the set of isometrics at the point from where the 
apparatus is calibrated. 

The indeterminacy would be lifted if either the Physics 
or the manifold itself provided a fundamental length. 
General Relativity does not have one built in; and space- 
time, being fundamentally featureless, cannot furnish one 
either. 

The method in this paper has both advantages and dis- 
advantages respect to the ones available in the literature. 
The most prominent difference is the use of coordinate- 
dependent quantities {e.g. Eq. to find the solution 
of a coordinate-independent equation (Vg = 0). Two 
rewards are worth of notice. First, the solution is gen- 
eral enough to cover the demands of Physics. Second, 
one can truncate the series to compute all of the metrics 



compatible with the geodesies up to the accuracy of the 
instruments. 

The results of this paper open the possibility of fur- 
ther studies on the stability of the inversion procedure. 
Such studies have been done for 3D Riemannian met- 
rics [3l| but there is no definitive answer for 4D pseudo- 
Riemannian ones. 
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APPENDIX 

Let us accept as hypothesis that an experimentalist 
has collected enough data to claim that the trajectories 
in a region of spacetime are given by 

x^' = {t,z) = (^t,v^t + ^at^ + ^(3t^ + 0{t'^)^ , (A.l) 

where Vz is constant in a trajectory, a and (3 arc constants 
for all trajectories, and the fitting is good to the order 
shown by O. For this parametrization, we get 

x^t=(^l,vz + at+^f]t^ + 0{t^)y (A.2) 
x^'tt = {0,a + pt + 0{t^)), (A.3) 
and Eq. ^ yields 

a + fSt + x^.x", (r^^,- 

(^Vz+at+^pf'^T\,,^ ^0{t'). (A.4) 

Since Vz is a parameter that names the geodesic, the 
Christoffel symbols should not depend on it. Thus, the 
coefficients of all powers of Vz must vanish in the last 
equation. This is equivalent to generate several equa- 
tions for different values of Vz ■ After solving the resulting 
equations, one gets 

^r\,^T\,^T\,^(, + 0{t'), (A.5a) 

^T\, = = T\, =6 + (t') , (A.5b) 
r\t = -a-l3t + 0{t^). (A.5c) 

The components not shown vanish. The unknown func- 
tions make the covariant vector anticipated by Levi- 
Civita's theorem. 
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The parametrization in Eq. (jA.ip implies that we 

should keep track of only up to O ((t, z)^). We can 
then write 

= 60 + Cut + ^i2Z + O {{t, zf) , (A.6) 

6 = 60 + 6it + 622 + O {{t, z f) , (A.7) 

where ^ij are constants. 

Let us accept as hypothesis that the rods and clocks 
used were such that if an element u of the tangent space 
at (0, 0) is acted by an element of the set S defined below, 
then the length of u remains constant. 

r / cosh(V') I sinh(V') V p re\ / A Q^ 
sinh(V') cosh(V^) (A-8) 

The boundary condition induced by S is 

50 = T Q 4 ) ' (^-^^ 
I 



where a will become the unknown constant conformal 
factor. 

We are set now to apply Ricardo's formula. The in- 
termediate calculations are straightforward yet cumber- 
some. The resulting metric does not always give back the 
connection of Eq. (jA.5P , but that was expected since not 
every value of yields a metric connection. It can be 
proven that the following constrains make the connection 
a metric one: 

62=6060, (A. 10) 

61=6060, (A.ll) 

61 = -"60 + 60 ^ • (A. 12) 

The components of the metrics compatible with the 
geodesies of Eq. (|A.1|1 are 



gtt = a- 



9 - 366oi - I860Z + - (8a2 + 54a6o - 1806'o " 816'o + 8I62) - 46o (2a -I- 276o) tz 

+ (46'o - 186'o - 962) z'] + O {{t, z)') , (A.13) 



gtz = a- 



(4a + 96o) t + 46oz + (-16a6o + 12^?o " 186'o " 862) tz - 2 (5a6o +13 + I86060) t 



166o6oz^ 



+ 0{it,zf) ,aiid (A.14) 



9zz = o- 



4 + 86o< + 166oz + (-8a6o + 12^?o - 862) + 486o6oi2 + 8 (46'o + $. 



62 z 



(A.15) 
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